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Abstract 

This paper considers the problem of selecting a subset of nodes in a two-hop wireless network to act 

as relays in aiding the communication between the source-destination pair. Optimal relay subset selection 

■ with the objective of maximizing the overall throughput is a difficult problem that depends on multiple 

H ■ 

' factors including node locations, queue lengths and power consumption. A partial decode-and-forward 

strategy is applied in this paper to improve the tractability of the relay selection problem and performance 
of the overall network. 

Note that the number of relays selected ultimately determines the performance of the network. This 
paper benchmarks this performance by determining the net diversity achieved using the relays selected 
and the partial decode-and-forward strategy. This framework is subsequently used to further transform 
relay selection into a simpler relay placement problem, and two proximity-based approximation algorithms 
are developed to determine the appropriate set of relays to be selected in the network. Other selection 
strategies such as random relay selection and a greedy algorithm that relies on channel state information 
are also presented. This paper concludes by showing that the proposed proximity-based relay selection 
strategies yield near-optimal expected rates for a small number of selected relays. 

Keywords - Greedy algorithms, partial decode-and-forward, superposition coding, relays. 
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1 Introduction 

Relay-assisted communication is a promising strategy for both centralized and decentralized communi- 
cation networks [1,2]. Two- hop relay-based communication is having a considerable influence on emerging 
standards both in local area networks, IEEE 802.11s [1] and broadband wireless access networks, IEEE 
802. 16j [2]. Two-hop relay systems consist of a source, a destination and one or more relays where the relay 
nodes work together as a single set of intermediaries between the source and the destination [3] . Direct trans- 
mission occurs between the source and the destination, and the relays assist the source only if the destination 
cannot decode the direct transmission. There are multiple concrete benefits of introducing these intermediate 
relays, which include improved system throughput and greater coverage [2]. Multihop relaying [4,5] is a key 
enabling technology for networks of the future, but before the performance tradeoffs of multihop relaying 
can be characterized, it is critical that the issues facing two-hop relaying be fully understood. 

Given that the source can enlist multiple nodes to simultaneously act as relays, two questions naturally 
arise. First, how many relays must the source enlist to aid its transmission to gain the maximum advantage 
for the resources consumed? Second, which of the nodes in the pre-existing network must be enlisted to 
act as relays? When multiple-relay selection is allowed, there are numerous tradeoffs that govern system 
performance [5-7] . While selecting a large number of relays offers the benefit of coherent combining, resulting 
in increased throughput and thus higher overall quality of service, it suffers from drawbacks as well. Firstly, 
system resources are drained faster when multiple relays are selected. Second, there are complexity and 
implementation issues - it is difficult to synchronize the transmissions from multiple disparate relays [8-11], 
and receiver complexity increases with the number of relays. A single relay can be selected to assist the source 
transmission [12-19], which offers lower gains in terms of total diversity and rate but is simpler to implement 
and consumes less power over the entire network. This paper has two goals. One goal is to understand 
the fundamental limits of multiple-relay selection to benchmark various relay selection algorithms. To this 
end, we focus on minimizing relay power consumption and treat implementation issues and complexity as a 
secondary concern. 

Regardless of the number of relays selected, it is difficult to determine which node(s) in the network must 
act as relays to aid the source transmission. For example, selecting the relay with the best channel to the 
destination may not be an optimal strategy, as this relay may be heavily loaded with traffic and running low 
on resources. Thus, relay selection is a very difficult problem, as selecting the "optimal" subset from the set of 
candidate relay nodes is affected by the presence of multiple parameters that govern system performance. In 
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particular, relay node selection often translates to a combinatorial optimization problem [31], which currently 
does not have an elegant polynomial-time algorithmic solution. 

The second goal of this paper is to provide algorithms for relay node selection that serve as a good 
approximation to the problem of optimal relay selection from the point of view of throughput maximization 
with power allocation. Moreover, we desire the algorithms to have low complexity and be highly intuitive in 
terms of design. Note that any selection algorithm is closely coupled with the transmission strategy employed 
in the network (such as decode/amplify/compress- and- forward). Thus, we discuss the transmission strategy 
employed in this paper and then delve into the details of the algorithms. 

In our paper, we use a partial decode- and- forward transmission strategy proposed in [2lj3. Partial 
decode-and-forward as described in [21] relies on a two-level superposition coding strategy introduced by 
T. Cover for broadcast channels [26] and further studied in [27-29]. Under this setting, the transmitter 
employs a layered coding strategy, allowing the receiver to decode the transmitter's message partially if it 
is incapable of determining it in its entirety. Note that the conventional decode-and-forward strategy as 
in [22] is a special case of the partial decode-and-forward strategy, and therefore, partial decode-and-forward 
is a useful tool that has all the properties of decode and forward incorporated into it. In particular, partial 
decode-and-forward offers both the diversity advantages of amplify-and-forward and the inherent robustness 
to noise of decode-and-forward [20]. 

The other main advantage of partial decode-and-forward is the tractability it lends to the relay selection 
problem. While multiple-relay selection based on partial decode-and-forward transmission does not read- 
ily lend itself to practical implementation, the resulting problem tractability facilitates the determination 
of valuable performance benchmarks, especially in terms of diversity gain. Our first contribution is the 
derivation of both the diversity gain and the generalized diversity gain that is achieved by allowing m of the 
candidate relays to assist the source. The resulting diversity analysis extends the single-relay result in [21] 
and further highlights the performance benefits of multiple-relay selection. We stress that the derived diver- 
sity gain is at most the diversity achieved by selecting m relays out of Kr candidate relays. For example, 
selecting the relay with the best end-to-end path between the source and the destination can yield a diversity 
gain of Kr + 1 [16]. 

We mention here that generalized diversity, which arises from the notion of generalized degrees of freedom 
^Note that this notion of partial decode-and-forward is distinct from the one in [22] as it is inspired by outage capacity. It is 
based on the superposition coding strategy for broadcast channels in [25], while the partial decode-and-forward strategy in [22] 
is derived from block-Markov coding. 
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[32,33], refers to the diversity achieved when the candidate relays have different transmit SNR values than 
that of the source. In our paper, we consider a specific case of generalized diversity where the SNR for each 
candidate relay is an exponential scaling of the source SNR, i.e. (Pj/cj^) = {Pt/cr'^)^ for relay i. This allows 
for performance benchmarking of networks where the relays may be operating on a different power budget 
than that of the source, including relay-assisted cellular networks. For example, the case where k < 1 can be 
modeled by a base station being assisted by mobile devices that are not currently handling their own voice 
traffic. On the other hand, the case where A; > 1 can be modeled by a battery-powered mobile device being 
assisted by fixed, dedicated relay nodes that are connected to a continuous power source. 

Our second contribution entails using the partial decode-and-forward framework as a platform to trans- 
form the relay selection problem into a relay-placement problem, whose solution suggests the "best" set of 
relays to be selected. There are two approximation steps: we first approximate selection of m relays by the 
problem of finding the m relays that are closest to a rate-maximizing location, and we show that obtaining the 
rate-maximizing location is equivalent to maximizing a signomial function [30] . Since signomial programs 
are, in general, not easy to solve, we further approximate the relay selection problem by the problem of 
finding the relays that are closest to the rate-maximizing location in a three-node line network. Obtaining 
this rate-maximizing location is equivalent to maximizing a polynomial over a given range of values, which 
can be accomplished using deterministic polynomial-time algorithms. The above polynomial approximation 
motivates two proximity-based algorithms (which we call Multiple Fan Out and Single Fan Out, detailed in 
Section [5|) that select relays based on their proximity to one of the rate-maximizing locations. In addition, 
we present a greedy selection algorithm (which we call Best Gains, also detailed in Section [5]) that chooses 
relays based on their channel gains to the destination and the amount of the source message that they have 
decoded. Here, the selected relays must have decoded at least one of the two messages from the source. 
Finally, we present a selection algorithm that randomly selects relay nodes (which we call Random Relays, 
also detailed in Section [5]) and compare the performance of all four algorithms - Multiple Fan Out, Single 
Fan Out, Best Gains and Random Relays. 

This paper is organized as follows. In Section II we describe the system model and introduce the two-level 
superposition coding strategy that will be used throughput the paper. The diversity analysis is shown in 
Section III. In Section IV, we present the analytical formulation of the relay selection problem and obtain a 
closed-form expression for the rate-maximizing relay position in a three-node line network. We present our 
proposed selection algorithms in Section V. We present simulation results in Section VI and conclude the 
paper in Section VII. 
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2 System Model 

First, we introduce the notation used throughout the paper. E denotes the mathematical expectation 
operator and In(-) represents the natural logarithm function, exp(-) represents the exponential function and 
r(-) is the gamma function. SNR represents the transmit-side signal-to-noise ratio at the source node. P{A) 
denotes the probability that an event A occurs. / (x) denotes the derivative of a function / with respect 
to its argument x. \\A\\ denotes the cardinality of a set A. |zp denotes the absolute square of a complex 
number z. f{x) ~ g{x) for large values of x represents the fact that f{x)/g{x) — > 1 as x ^ oo [20, Pg. 3067]. 

Consider the two-hop wireless network in Fig. [TJ The network consists of a single source t, a single 
destination r and Kr relays interspersed throughout the region between the source and the destination. Let 
di^n denote the distance between nodes i and n. Let hi^n denote the channel between nodes i and n. 

2.1 Key Assumptions 

We make the following critical assumptions in this paper: 

• Each relay operates in a half-duplex mode and employs a single antenna. 

• Additive white circularly symmetric complex Gaussian noise nj with mean and variance o"^ is present 
at each receiving node i during time slot k. 

• \hi^n\ is a Rayleigh-distributed random variable. Thus, the real and imaginary components of h^ i^ are 
mutually independent Gaussian-distributed random variables, each with mean and variance (1/2) • 
E(|/ij_„|^). This assumption simplifies our analysis and is typically used in the literature to obtain 
insights on the performance of real-world wireless systems. 

• Our transmission strategy has arbitrarily long codewords that are generated using an i.i.d. Gaussian 
distribution of suitable variance to meet the overall power constraint. Note that the capacity and 
thus, the capacity achieving input distribution (if any) of the additive Gaussian noise relay channel is 
in general unknown. Thus, we choose this coding strategy for two reasons: 1) it has been found to 
be optimal in most of the special cases whose capacity is known (physically/reversely degraded [22], 
orthogonal channels [23] and uniform phase fading [24]) and 2) it yields explicit rate expressions and 
intuitive coding strategies. The destination and all potential relay nodes employ typical set decoding 
as defined in [34]. 
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• The source knows the exact channel state for ah of the channels in the network in Fig. [TJ Each relay 
knows the exact state of its channel from the source. The destination knows the exact state of its 
channel from each of the relays and from the source. 

• Time is slotted and the channel is constant in every time slot. 

• Each time slot is large enough to admit an arbitrarily small probability of error as long as the rate in 
that state is below the maximum achievable for that state (this is also referred to as the block fading 
assumption). 

• A log-distance path loss model is applied [35]. Let Ac, do, and fi denote the carrier wavelength, the 
reference distance, and the path loss exponent. Then, the channel gain between nodes i and n is 



2.2 Partial Decode-and- Forward 

All relays perform partial decode-and-forward operations based on the two-level superposition coding 
strategy in [21]. The source transmits xt^i during the first time slot, where 



and the source allocates power (3Pt to xi and power (3Pt to X2, where (3 G [0, 1] and P = 1 — (3. Note that 
xi and X2 are codewords from codebooks with elements that are generated i.i.d. according to zero-mean 
Gaussian distributions with variance pPt and pPt, respectively. 

The destination and all candidate relay nodes employ typical set decoding to decode xi and X2- The 
candidate relays and the destination initially attempt to decode xi. If node i can decode xi then it attempts 
to decode X2- Two channel thresholds, \hi\ and |/i2|, are chosen to determine the set of received rates for 
this two-level coding strategy. Then, xi can be decoded at the rate Ri [21], where 



(Ae/47rdo)2(d.,„/do)-^. 



(1) 



Xt,l = XI+X2 



(2) 




(3) 



while X2 can be decoded at the rate R2 [21], where 




(4) 
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Note that if node i attempts to decode xi or X2 at a higher rate than Ri or R2, respectively, the resulting 
probability of error is bounded away from zero. 

The received signals at the candidate relay i and at the destination during time slot 1 are, respectively 

yi,i = ht^iXt,! + (5) 

yr,l = ht^rXt,! + Ur^l. (6) 

If the destination can decode both xi and X2, it broadcasts this information to the entire network and the 
source prepares to send xt^2 during time slot 2. If the destination can only decode xi, or if it cannot decode 
xi, it broadcasts this information to the entire network. The source then selects a subset of the candidate 
relays to assist its transmission. 

For relay i, if \ht^i\ < then it cannot decode xi and it does not transmit during time slot 2. If 

l^il ^ \^t,i\ < 1^2!) then a selected relay i can only decode xi and will forward xi to the destination during 
time slot 2. If \ht^i\ > I/12I1 then a selected relay i can decode xt,i and will forward xt^i to the destination. 
Thus, relay i allocates power Pi to its transmission Xr,i to the destination, where [21] 

if \ht^i\ < \hi\ 

y^xi if 1^1 1 < \ht^i\ < |/i2[ (''') 

T^i^i + X2) if \ht,i\ > \h2\ 

For each relay i, we set [3i = (3 for the majority of this paper; in Section [6.11 we investigate the performance 
impact of varying /3j with respect to (3. 

Thus, if A denotes the set of all relays that transmit during time slot 2, the destination receives 

?/r,2 = ^ hi^rXr,i + ^^r,2 (8) 

during time slot 2. After time slot 2, if the destination can decode x^^i, the received rate is Ri + i?2- If the 
destination can only decode xi, the received rate is Ri, and if the destination cannot decode xi, the received 
rate is 0. Note that this two-level coding strategy can be generalized to a multiple-level approach based on 
broadcast strategies introduced for the single user and MAC channels [25]. Once the two-level strategies 
and algorithms are understood, their generalization to n > 2 levels is relatively straightforward but leads 
to extremely unwieldy expressions. Moreover, it is unclear if using a multiple- level approach will provide 
significant gains in performance. Thus, we have chosen to limit ourselves to a two-level transmission strategy 
in this paper. 
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Our proposed multiple-relay selection algorithms choose A to maximize the expected rate subject to a 
sum power constraint over all relays i A. In the next section we derive both the diversity gain and the 
generalized diversity gain via selecting m relays. 

3 Diversity Performance 

After m relays are selected to transmit to the destination during time slot 2, we consider the resulting 
diversity gains Ki{m) and K2{m) for Ri and i?2j respectively. Let Pout{Ri-,A) denote the probability that the 
destination cannot decode xi after time slot 2, and let Pout[R2-: A) denote the probability that the destination 
cannot decode X2 after time slot 2. 

For the diversity analysis, we set the relay powers Pi = Pt for i G {1, 2, . . . , m}, and so the SNR is Ptja^ . 
The diversity gains are obtained by observing that the outage probabilities Pout{P\^A) and PoutiP2-,A) are 
proportional to SNR^'^'^^'^^ and SNR~'^^^"^\ respectively as the SNR Pt/a'^ approaches infinity. We reiterate 
that these diversity gains are at most the diversity achieved by selecting m relays out of Kr candidate relays. 

Theorem 1. Selecting m relays to transmit during time slot 2 yields a diversity gain of 

Ki{m) = m + 1 

and K2{m) = Ki{m). 

Proof. See Appendix |X1 □ 

We note that obtaining the diversity gain of m + 1 entails a relatively straightforward extension of the 
single-relay analysis for decode-and- forward relaying in [20, Section IV. B]. In the two-level transmission 
strategy that we consider, the destination still attempts to decode the transmission from the source even if 
at least one relay fails to decode xi or X2- Note that the single-relay analysis in [20, Section IV. B] ignores 
the direct link transmission if the relay makes a decoding error, and so a direct extension of the analysis 
in [20, Section IV. B] would yield a diversity gain of m instead. 

We also perform a generalized diversity analysis where the relay powers are such that (Pj/cr^) = (Pt/a^)'' 
for i £ {1, 2, . . . , m}, where A; is a real number. The generalized diversity gains and Kfl^T-) are obtained 

by observing that the outage probabilities Pout{Ri, A) and PoutiR2, A) are proportional to SNR^'^i^'^^ and 
SNR~'^^^^\ respectively as the transmit-side SNR Pt/a"^ approaches infinity. 
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Theorem 2. Selecting m relays to transmit during time slot 2 yields a generalized diversity gain of 

Kf(m) = 



km + 1 \i k <1 
m + 1 if > 1 



and ^2(711.) = K2(?n.). 

Proof. See Appendix [Bj □ 

The generalized diversity gain K^{m) = ^2(712) has the following intuitive interpretation. If /c < 1, each 
relay is no better than the source in terms of transmit power, so the worst-case error event is determined 
by all of the relay-to-destination channels. In particular, this event occurs when all m suboptimal relays 
attempt to forward xi or X2 to the destination. On the other hand, if A; > 1, each relay is better than the 
source in terms of transmit power, so the worst-case error event is determined by all of the source-to-relay 
channels. In particular, this event occurs when all m superior relays cannot decode either xi or X2- 

4 Rate-Maximizing Relay Position 

We formulate the relay selection problem for an arbitrary number of selected relays, and then we show 
how this problem can be simplified by considering a three-node line network. 

4.1 Optimal Relay Placement in General Network 

Consider the case where a subset A of the available relay nodes {1,2, Kr} are selected to assist the 
source. Let h denote the channel between a transmitting node and a receiving node. The received rate at a 
receiving node via decoding xi is [21] 

and the received rate at a receiving node via decoding X2 after decoding xi is [21] 



cm')^in{i + ^-!fi^). (10) 



The expected rate of the two-level superposition coding strategy is 

RscA^) = (1 - PoutiRi,A))Ri + (1 - Pout{Ri,A))il - Pout{R2,A))R2 (11) 
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and so the relay selection problem can be formulated as follows 

max Rsc2iA) (12) 

AC{l,2,...,Kr} 

subject to ^ Pj < Pmax and < Pj < Pi^max Vi E A. 

It is apparent from (jl2p that the relay selection problem is also a power allocation problem. In particular, 
if a relay i is not selected, it is assigned a power Pj = 0. On the other hand, if a relay i is selected, it is 
assigned a power Pj > according to the solution to (fT^ . 

Let A denote the set of all relays that cannot decode xi, and let Q denote the set of all relays that can 
decode xi but cannot decode X2- The probability that the destination cannot decode xi after time slot 2 
can be obtained by generalizing [21, (13)] as 

Pout{RuA) = E ( E [[X{P{Ci{\K5?)<Ri)\ 



X 



\{P{Ci{\ht,e?)>Ri,C2{\htM?)<R2)\[ n P{C2{\ht,r,?)>R2) 



Each term in the inner sum in (|13|) represents a scenario where a selected relays cannot decode xi, ^ selected 
relays can decode x\ but cannot decode X2, and the remaining ||^|| — a — ^ selected relays can decode x^- 

The expressions in (|13p are fairly involved, so we consider the high-SNR regime for ease of analysis. In 
Appendix El we prove that 

P(Ci(|Vp)<PO ^ (14) 

•-^t.^ (^t/f^ ) X (1 - /3exp(Pi)) 

P(Ci(|/it,e|') >Pi,C2(|/it,e|') <P2) ~ 1 (15) 
P(C2(| VI') > P2) ~ 1 (16) 



and 



/ A IVP/^^f + E^g(^\(Aue)) |/ie,,pPe \ ' 

1^ 1^ l^*..|W + E,e(^\(Aue))IVP/?^. + a2+Z^ ^2 ^ ^ 



\ ll-AII-Q+i 
exp(Pi) — 1 ^ 

X 



(Pt/a2) X (1 - /?exp(Pi)) ; (Pli - ^ ^ ^^^^^^^ 



1 1 -|-|- 1 
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The probability that the destination cannot decode X2 after time slot 2 is 

Pout{R2,A) = E ( E ((llP{C2i\ht,s\^)<R2 

0<a<\\A\\ \ACA,\\A\\=a WseA 



X 



X 



n P(.C2(.\ht,ef) > R2)] 
.ee{A\A) / 

p(^C2(^\h,r\'+ E IVP^<^2^^y (18) 

Each term in the inner sum in (jlSp represents a scenario where a selected relays cannot decode X2 and the 
remaining ||^|| — a selected relays can decode X2- 

The expressions in (jlSp are also fairly involved, so we again consider the high-SNR regime for ease of 
analysis. In Appendix [Aj we prove that 



P{C2{\h,M') < R2) ~ 7^>< ~TTT-^ (19) 



1_ exp{R2) - 1 

PiC2{\ht,e\^) > R2) ~ 1 (20) 



and 



P[C2l\ht,r\^+ E IVrf) <^2] 



exp(i?2) - 1 \ 

m/'^') I (Pll - a + 1)! - Gl ^^f-J^^^ {Pe/Pt)Gl^ 



-a + Dl'' G?.^U.(PJP)Gl- ^^^^ 



It is apparent that (jl2p is an optimization problem with linear inequality constraints. Also, from in- 
specting (fT3]) - (fT7|) it is clear that Pout{Rii is a nonlinear function of Pj V i € ^ in the high-SNR regime. 
In addition, from inspecting (fT8|) - ([2T]) it is clear that Pout{R2, A) is a nonlinear function of Pj V i € ^ 
in the high-SNR regime. Then, the preceding analysis shows that in the high-SNR regime, Psc,2(-4) is a 
nonlinear function of Pj for i £ A. Thus, nonlinear programming techniques such as sequential quadratic 
programming [36] can be applied to solve (|12p in the high-SNR regime. 

The relay selection problem (jl2p can also be approximated as a relay placement problem where m relays 
in Fig. [T]are chosen to assist the source. The key idea behind the relay placement problem is to hypothetically 
place m relays in the locations that would maximize Psc,2("4). Then, the m relays in Fig. [T]that are closest 
to the rate-maximizing locations are selected. It is also assumed that each selected relay i employs the same 
power Pi = Pmax/m. 
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To solve for the rate-maximizing locations, recall from ([T]) that G?^ = (Ac/47r(io)^((ii,n/'^o) ^ = (c^i,n) ^X- 
Without loss of generality, assume that the source is located at (0,0) and the destination is located at (dt.r,0). 



If relay i is located at {ai,bi), then dt.i = y^f + bf and di^r = \J {dt,r — o-iY + bf. The rate-maximizing 
locations are 

{al,b\,...,al^,b*^} = arg max Rsc,2{A) 

ai,bi,...,am,,bm 



subject to ll^ll = m, ^ Pi < Pmax and < < Pi^max Vi G A. 



In particular, by considering the binomial series X^^Q(a + b)'' where k is a real number, we see that 



oo 



r(^/2 + 1) ^2k^,-2k ^22) 



^^A:!r(^/2 + l-A;) 



and 

- | .,r^2^-.) '*-- -°-)"r" (23) 

We assume that < Oj < dt^r for each relay i since the relays are interspersed throughout the region 
between the source and the destination. Also, assume without loss of generality that 6^ > for each relay i 
since dt^i and di^r are functions of 6?. Let the m selected relays be located at (oi, bi), (a2, 62), • • • > {d-m, bm)- 
Recah from ([13]), (HH), ([IT]), ([H]), and ([H]) that in the high-SNR regime, Rsc,2{A) is a function of 
^i~n ~ {di,n)^/x- Then, from ([22]) and ([23]), we see that Rsc,2{A) is a function of {oi, 61, ... , am, bm}- Since 
we have assumed that Oj and 6j are positive for each relay i, and the binomial coefficients in (j22p and (|23p are 
not necessarily positive, Rsc,2{A) is a signomial function [30] of {ai, 61, ... , a^, ^m} in the high-SNR regime. 

Signomial programs usually do not admit efficient solutions via geometric programming unless the objec- 
tive function and the associated inequality and equality constraints satisfy certain conditions [30]. Next, we 
show that given a three-node line network, Rsc,2iA) reduces to a polynomial function of the relay position d. 

4.2 Optimal Relay Placement in Line Network 

We consider a line network with Kr = 1. The source is located at (0,0), the destination is located at 
{dt^r,0) and the relay is located at (d, 0) where < d < dt^r- The outage probability PoutiRi, A) can be 
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written as [21] 

Pout{Ri,A) = P{Ci{\ht,i\^) < ) < 

+ P{Ci{\ht,i\^) > Ri,C2{\Ki\^) < R2) 

xP In 1 + - — -f^ + — '— < Ri 

+ PiC2{\ht,i\^) > i?2)p(ci(|Vf + '^''^'^' ) < (24) 
and the outage probability Pout (-^2! -4) can be written as [21] 

Pout{R2,A) = P{C2{\ht,l\^) < R2)P{C2{\ht,r\^) < R2) 

+ P(C2(|/ii,i|2)>i?2)p(c72(|Vf + %r^) <^2). (25) 

As in Section 14.11 the expressions in ()24p and (j25p are fairly involved, so we again consider the high-SNR 
regime for ease of analysis. 

In Appendix 1X1 we prove that (j24p simplifies to 

P (R A\ ( ^ exp(Pi)-l \( 1 exp(Pi)-l \ 

''^ - I ^2 (p^/^2) ^ (1 - ^exp(Pi)) ] \Gl (Pja^) X (1 - ^exp(Pi)) ) 



1 exp(Pi) - 1 \ / 1 exp(Pi) - 1 

iPi/Pt)Gl, {Pt/a^) X (l-^exp(Pi)) ) I (P^/aS) x (l-^exp(Pi)) 



and we prove that ([25]) simplifies to 



„ I exp Pa -1\ / 1 exp P2 -1\ 

+ 2"l^(Pi/p,)Gf/ m/-^) ) [gi mi<y') )' ^ 

Let Gi = (exp(Pi) - l)/((Pt/cj2) x (1 - ;3exp(Pi))) and G2 = (exp(P2) - l)/(,9(Pt/cT2)). Then 

RscM) = il-Pout{Ri,A))Ri + {l-PoutiRi,A)){l-Pout{R2,A))R2 

~ Pi X (1 - Gjx^ X + (P^/Pi) X (dt,, - d)^)) 

+ P2 X (1 - Gjx^ X + (P^/Pi) X (dt,, - dr)) 

X (1 - G^x' X + (1/2) X (Pi/Pi) X idt,r - dr))- (28) 

For integral values of the path loss exponent //, finding the rate-maximizing relay position d is equivalent to 
maximizing a polynomial over < d < dt^r- For example, if /i = 2, Psc,2(-^) is a fourth-degree polynomial in 
d. Maximizing Psc,2(-^) with respect to d is then equivalent to finding the roots of a cubic equation that lie 
in < d < dt^r, assuming that at least one exists. 
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5 Relay Subset Selection Algorithms 

The analysis in Section [4.11 shows that the combinatorial optimization problem (jl2p can be approximated 
by considering a given value ofmE{l,2,..., Kr} and maximizing a signomial function of the relay locations 
(oi, . . . , (om, bm) in the high-SNR regime, which yields the rate- maximizing set (oi, 6i), . . . , {am,bm)- 

Note that since Rsc,2{^) is a signomial function in the high-SNR regime, relays that are located close 
to any of the points in the rate-maximizing set should still yield high expected rates due to the inherent 
smoothness of signomial functions. This motivates the following proximity-based algorithm for solving ()12p . 

Algorithm 1. Multiple Fan Out 

Step 1: For a given value of m E {1,2, Kr}, maximize -Rsc,2(-4) over all relay locations 
(oi, . . . , {am, bm) to find the rate-maximizing set (ai, . . . , {am,bm)- 

Step 2: Set i = 1 and A = (l). 

Step 3: For relay n, where 1 < n < K^, compute d{n) where d{n) is the distance from relay n to 
{ai,bi). If relay n is at location {an,bn), d{n) = \J (an — ai)^ + {bn — biY- 

Step 4: Find the closest relay n to {en, hi) not in A and let ^ = ^ U {n} and P„ = Pmax/^- 

Step 5: If ll^ll = m, stop. Otherwise, let z = i -|- 1 and return to Step 3. 

We call the above relay selection algorithm Multiple Fan Out because the process of relay selection is 
analogous to a search party fanning out from its initial location. Here, the objective is to "fan out" from 
(ai, bi), . . . , {am, bm) until m relays have been selected. 

Note that in Step 2, we set i = 1 and increase i in Step 5. It turns out that ai = ■ ■ ■ = am and 
5i = • • • = 5m in Step 1, so the initial assignment of i in Step 2 and its iteration in Step 5 are irrelevant. 
The Multiple Fan Out algorithm, then, reduces to finding the m closest relays to the single rate-maximizing 
point (ai, 

Step 1 involves maximizing a signomial function, which usually does not admit an efficient solution. To 
obtain a more tractable problem, the analysis in Section 14.21 shows that in the case of a three-node line 
network, Rsc,2{A) is a polynomial function of the relay location d in the high-SNR regime. Maximizing 
Rsc,2{A) yields the rate-maximizing relay location d. 

Also, since Rsc,2{A) is a polynomial function in the high-SNR regime for a three- node line network, relays 
that are located close to the rate-maximizing relay location d should still yield high expected rates due to the 
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inherent smoothness of polynomial functions. This motivates another proximity-based algorithm for solving 
(jl2p . Again we assume that the path loss exponent /x takes on an integral value. We also assume that exactly 
m relays are to be selected, which further simplifies the algorithm. 

Algorithm 2. Single Fan Out 

Step 1: Maximize (j28p to find the rate-maximizing relay location (d, 0). 

Step 2: For relay n, where 1 < n < Kr, compute d{n) where d{n) is the distance from relay n to 
(J, 0). If relay n is at location {an,bn), d{n) = [an — dY + ^n- 

Step 3: Sort the set of relays {1,2,..., Kj.} as {ai, a2, . . . , a^^}, where 

d{ai) < d{a2) < ■ ■ ■ < d{aKr)- 

Step 4: Find the closest relay n to (d, 0) not in A and let ^ = ^ U {n} and P„, = Pmax/'>^- 
Step 5: If ll^ll = m, stop. Otherwise, return to Step 4. 

We call the above relay selection algorithm Single Fan Out because d is computed via analysis of a single- 
relay line network. Note that both the Multiple Fan Out and Single Fan Out algorithms are greedy strategies 
in that the order of procession through the list of relays is based on their proximity to (aii, yi), . . . , {x^^ Vm) 
and d, respectively. Greedy algorithms are useful for the problem at hand in that they possess an inherent 
simplicity, and their run-times are usually simple to characterize. 

We propose another greedy approach for selecting A. For simplicity, we assume that at most m relays 
are to be selected. 

Algorithm 3. Best Gains 

Step 1: For relay i, where 1 < i < K^-, compute |/ii,rP- 

Step 2: Sort the set of relays {1, 2, ... , Kj.} as {oi, 02, ... , axr}-, where 

Step 3: Let i = 1 and ^ = 0. 

Step 4: If relay Oj has decoded xi, then ^ = {aj} and Pj = Pmax/fn. 

Step 5: If ll^ll = m or i = K^, go to Step 6. Otherwise, let i = i -|- 1 and return to Step 4. 

Step 6: If ll^ll < m, let Pi = Pmax/\\A\\ for each relay Oj G A. 
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We call the above relay selection algorithm Best Gains because the order of procession through the list 
of relays is based on their channel gains to the destination. The objective is to choose relays that will be able 
to reliably transmit to the destination during time slot 2. Note that a check is performed on each selected 
relay in Step 4 to ensure that it will be able to forward at least xi to the destination. 

To obtain a lower bound on the performance of the above greedy algorithms, we propose the following 
algorithm whereby relays are randomly selected to transmit during time slot 2. 

Algorithm 4. Random Relays 
Step 1: Let ^ = 0. 

Step 2: Randomly select a relay i € {1, 2, ... , Kr}\A and let A = AU{i} along with Pi = Pmax/m. 
Step 3: If ll^ll = m, stop. Otherwise, return to Step 2. 

Since the destination employs a diversity combining approach to receive the signals from all of the selected 
relays, the performance of the Random Relays algorithm should approach that of the other proposed relay 
selection algorithms as m increases. 

6 Simulation Results 

6.1 Performance of Relay Selection Algorithms 

We place the source at (0,0) and the destination at (100,0). We use the Worldwide Interoperability for 
Microwave Access (WiMAX) signaling bandwidth of 9 MHz [38], and given a noise floor of -174dBm/IIz this 
yields a noise value = — 104dBm. We also have a carrier frequency fc = 2.4GIIz along with a reference 
distance do = Im and a path loss exponent ^ = 3. We randomly place Kr = 20 relays in the region between 
the source and the destination. 

Fig. [2] shows how the expected rate Rsc,2iA) varies with the number of selected relays ||^|| = m for the 
algorithms that we have proposed. Here we fix the source's power Pt = 6dBm and the relay sum power 
constraint Pmax = Pt- The fraction of the source's power allocated to xi is /? = 0.75 and the decoding 
thresholds for xi and X2 are \hi\ = 7.4 • 10~^^ and |/i2| = 1.25 • lO"^'^, respectively. In the case of the Best 
Gains algorithm, we only consider cases where the number of selected relays |[^|| = m. We obtain the 
rate-maximizing set (oi, 6i), . . . , {am,bm) for the Multiple Fan Out algorithm via the fmincon function from 
Matlab, which employs a sequential quadratic programming method. 
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We see that the greedy Best Gains algorithm yields the highest expected rate for all of the proposed 
selection strategies. This is due to the fact that the Best Gains algorithm biases relay selection towards 
those relays that have good channel gains to the destination and can also transmit during time slot 2; 
this minimizes the chances of an outage event occurring at the destination where it cannot decode xi. 
On the other hand, the Fan Out algorithms select relays that are close to ergodic rate-maximizing points 
without considering their decoding status and their instantaneous channel gains to the destination. Thus, the 
Best Gains algorithm attempts to optimize relay selection for each source transmission, though additional 
overhead is incurred relative to the Fan Out algorithms since the relays must inform the source of their 
decoding status and their channels to the destination. 

Also, the Single Fan Out algorithm offers virtually the same performance as the Multiple Fan Out al- 
gorithm, which demonstrates the utility of our simplifications of the relay selection problem. Here, the 
relays that are close to the rate-maximizing set for the Multiple Fan Out algorithm are also close to the 
rate-maximizing position (d, 0) for the Single Fan Out algorithm. In addition, as the number of selected 
relays ||^|| increases, each strategy yields a higher expected rate which approaches the maximum expected 
rate. Finally, note that the performance gap between all of the proposed strategies decreases as the num- 
ber of selected relays increases. This is due to the fact that selecting multiple relays yields an SNR gain 
at the destination that gradually overcomes the loss from selecting relays that might not be close to the 
rate-maximizing positions that are computed by the Multiple Fan Out and Single Fan Out algorithms. 

Fig. [3] shows how the expected rate -R<jc,2(-^) varies with the number of selected relays ||^|| = m for two 
relay power allocation strategies. We use the same system parameters as in Fig. [2l except that we randomly 
place m relays in the region between the source and the destination instead of Kr = 20 relays. The Optimal 
Power Allocation strategy entails solving the relay selection problem in (jl2p , and the Equal Power Allocation 
strategy assigns equal power to all of the selected relays. We also set Pmax = Pt- 

We observe that the Equal Power Allocation strategy offers comparable performance to the Optimal Power 
Allocation strategy. This illustrates the utility of low-complexity strategies that reduce the computation time 
inherent to interior-point methods that are needed to solve the optimization problem (|12p . 

Fig. m shows how the expected rate Rsc,2i-^) varies with the average received SNR at the destination 
for different ratios between the relays' and source's powers. When the average received SNR values at the 
destination are OdB, 2dB and 4dB, the source's power takes on values Pt = — 6dBm, Pt = — 4dBm and 
Pt = — 2dBm, respectively. 

We see that as the average received SNR at the destination increases, the expected rate increases for each 
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value of Pi I Pi. Note that for a fixed value of Pi/Pt, increasing the average received SNR at the destination 
entails increasing Pi and Pt- For a fixed value of the average received SNR at the destination, the expected 
rate decreases as Pi/Pt decreases, which corresponds to a decrease in Pj. Thus, even though the three selected 
relays yield an SNR gain at the destination in time slot 2, this gain decreases as the relays' power decreases. 

Fig. [5] shows how the expected rate Rsc,2{-^ varies with the average received SNR at the destination for 
different values of the relays' power split We set (3 = 0.75. 

We see that as the relays' power split /3i decreases, the expected rate increases for all average received 
SNR values at the destination. Note that as f3i decreases, /3j increases, which leads to an increase in R2 as 
seen in On the other hand, as f3i increases, ^ shows that Ri decreases. Fig. [5] shows that the increase 
in i?2 overcomes the decrease in Ri. 

7 Conclusion 

We have studied the problem of selecting a set of relay nodes to forward data in a two-hop wireless 
network. We have considered a scenario where all relay nodes perform partial decode-and-forward operations 
based on a superposition coding strategy. For this setup, we have shown that relay selection can be initially 
approximated by the problem of finding the relays that are close to a rate-maximizing location. Finding 
the rate-maximizing location is usually computationally intensive, so we further simplify the relay selection 
problem by solving for the rate-maximizing location in a three-node line network. These results motivate two 
proximity-based relay selection algorithms, where relays are chosen to forward data based on their proximity to 
one of the rate-maximizing locations. We also demonstrated that the proximity-based algorithms outperform 
a random relay selection algorithm and yield rates close to those yielded by a greedy strategy that is based 
on channel state information. In addition, we derived the diversity gain achieved by having multiple relays 
assist the source. We also illustrated the performance impact of varying system parameters such as the ratio 
between the relays' and source's powers. 

As noted in the Introduction, selecting the optimal subset of candidate relay nodes to assist a source 
is a difficult problem, and the proposed selection strategies are mainly intended to offer key insights. In 
particular, the proximity-based algorithms motivate intelligent relay placement in a general two-hop static 
network with non-Rayleigh fading. System designers can experiment with different network topologies and 
determine a throughput-maximizing configuration, where the achieved throughput would depend on the level 
of interference between the transmissions from distinct relays. Also, the information-theoretic analysis in this 
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paper can be modified to support more practical transmission strategies. By applying cutting-edge coding 
strategies such as punctured low-density parity-check (LDPC) and turbo codes, the superposition coding 
approach that is employed in this paper can form the basis of a hybrid- ARQ strategy in a multihop network. 



A Proof of Theorem [T] 

The probability that the destination cannot decode xi after time slot 2 is 
Pout {Ri , ■A) = 



E ( E (llP(Ci{\h,s\')<Ri)) 

■.,^<m),a+^<m V ACyl.ec^, || A||=a,||e||=5, A f] 6=0 V <5gA / 

llP{Ci{\ht^e\')>Ri,C2{M)<R2)]( n P{C2{\ht,r,\^) > R2)]] 
eee / \nGU\(Aije)) / / 



,,?G(^\(Aue)) 



Each term in the inner sum in (129^ represents a scenario where a selected relays cannot decode xi, ^ selected 
relays can decode xi but cannot decode X2, and the remaining m — a — ^ selected relays can decode X2- 
Note that for a Rayleigh fading channel h, 

P{C,{\h\') < Ri) = pflnfl+ ly/^* ^ <R, 



(30) 



\h\^pPt + a^ 

>{Ri) - 1 
P exp(i?i) 
exp(iii) — 1 



Pl |/i|2 < exp(i?i) - 1 ^ 
' l-/3exp(i?i) 



E(|/i[2) (l-/3exp(i?i))Pt/a2 
where ([30]) follows from [20, Fact 1]. 

Also, for a Rayleigh fading channel h, 

~ 1. (31) 

In addition, for independent Rayleigh fading channels hi and /12, note that 

= P(Ci(|/li|2 + |/i2[2) <Pi) 

= p(\h,f^\h2?< -MRi)-^ 

V l-/3exp(Pi) Pt, 

exp(Pi)-l Y ^32) 



2E(|/ii|2)E(|/i2|2) V(l - ;9exp(Pi))Pt/a2 
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where ([32]) follows from [20, Fact 2]. 

Also, for a Rayleigh fading channel h, 

P{C2{\h\^) > R2) - 1. 



In addition, for independent Rayleigh fading channels hi, /12 and /is, note that 



+ \h2\^l3Pt + \hs\'^l3Pt+a2 
P(Ci(|/li|2 + |/l2|2 + |/l3|2)<i?i) 

1 / exp(Ri)-l 

6E(l/iil^)E(l/i2P)E(lfe3l^) V(l-/3exp(i?i))Pt/<72 



where p4p follows from [6, Appendix B]. 



We use ([SOD, ((MI, (l32|), ([Ml) and dMl) to see that 

I2/3D I I/, |2 



- 1^ 1^ \Kr?m + Y..^^A\A)\Kr?~PPt + <y^) ' 



P(cil\ht,r\^+ K,rn<Rl] 

\ \ ue(A\A) / / 



e{A\A) 

1 / exp(i?i)-l \-(— +1) 1 



n 



(m - a + 1)! E(| VP) V(l " /?exp(Pi))Pi/aV ^^11^^ E([ VP) 

Thus, the high-SNR behavior of Pout{RiT^) is 



and so we obtain a diversity gain of Ki(m) = m + 1 for decoding xi at the destination. 
The probability that the destination cannot decode X2 after time slot 2 is 



Pout{R2,A) = E ( E ( n^(C'2(|V|2)<P2)) 

0<Q:<m \ ||A||=a,ACyt \<5eA / 



X 



n P{C2{\ht,e\^) > R2) 
~0e{A\A) / 



xp[c2[ivi'+ E ivn <^2) J 

V V eeU\A) / / / 



ee(^\A) 
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Each term in the inner sum in (1370 represents a decoding scenario where a selected relays cannot decode X2 
and the remaining m — a selected relays can decode X2- 
Note that for a Rayleigh fading channel h, 



PiC^m') < R2) = P(^ln(^l + ^^) <i?2) 



p(|^|2^exp(i?,)-l^a^ 

1 ^ eMR2) - 1 ^3g^ 



E(|/i|2) ^P,/a2 
where ([MD follows from [20, Fact 1]. 

Also, for independent Rayleigh fading channels hi and h2, note that 

P{C2{\h,\^ + \h2\')<R2) = ^(ln(l + ^^ + ^^)<i?2 



/3 

1 / exp(fi2) - 1 

2E([/li|2)E([/l2|2) V /3Pi/fT2 

where follows from [20, Fact 2]. 

In addition, for independent Rayleigh fading channels hi, /12 and h^, note that 

V fi Pt 

3 



where (|40p follows from [6, Appendix B]. 

We use das]), (PH]), ([39]) and (00]) to see that 



X 



6'e(^\A) 
/exp(i?2)-l^"^™~°+'^ 



n 



(m-a + l)! E(IVP) V m/<r' J .gf/vA) ^^'^-'^'l') 

Thus, the high-SNR behavior of Pout{R2, is 

PoutiR2,A) ~ 1 (4 x(ir-"x 4 



-(m+1) 



C72 



(39) 



6E(|/ii|2)E(|/l2|2)E(|/l3l2)V J ^ ' 



p(c2[|Vl'+ \he,r?]<R2] (41) 



(42) 
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and so we obtain a diversity gain of K2{m) = m + \ for decoding X2 at the destination. 

Thus, we conclude that selecting m relays allows us to reap a diversity gain of m + 1 for both Ri and i?2- 



B Proof of Theorem [2] 

The proof of Theorem [2] is similar to that of Theorem [1] in Appendix 1X1 

First, we consider the decoding of xi at the destination. Recalling that {Pi/a"^) = [Pt/cr'^)^ for each relay 
i, we can use (p5]) to see that 

- [ [ \ht,r\'m + EueiA\A)K,r\'PPu + <^') V 



1 1 / exp(iii) - 1 



X 



-(m—a+l) 



n 



(m - a + 1)! E((Pi/a2)| V|2) V(l - /5exp(i?i))y ^^Jl^^ Ei{Pj a^)\K^r\^) 

1 1 / exp(iti) — Pt \ TT 



n 



(m-a + l)! E(lVrP) V(l - ^exp(iii))y \^a2y ueiA^)^^^^"''^^^' 

Thus, for a given integer value of a € {0, ... , m}, the high-SNR behavior of Pout{Ri, ■A) is 

~ ((-|j j x(l)'3x(ir-"-'3x 

(\ -(fcm+l+«(l-fc)) 

We then minimize km + 1 + a{l — k) over all a G {0, . . . , m} to obtain the generalized diversity gain K^{m) 
in Theorem [2j 

We then consider the decoding of X2 at the destination. Recalling that {Pi/a"^) = {Pt/cr'^)^ for each relay 
i, we can use (HTl) to see that 



eg(yt\A) 



1 1 / exp(fi2)-l V^"'"+'^ T-r 1 

_ a + 1)! E{{Pt/a^)\ht,r\^) \ J }■} E{{P, / a^)\K,r\^) 



1 If exp(i?2) - 1 \ ( -f* \ 



n 



(m - a + 1)! " E(| VP) 'A P J W) ^J^}^^^ ^(1^12^ ' 
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Thus, for a given integer value of q € {0, . . . ,m}, the high-SNR behavior of Pout{R2,^A.) is 

(/ \ -1\ a / \ —k{m-a)+l) 

i^l ] x(ir--x f^j 

We then minimize km + 1 + a(l — k) over all a G {0, . . . , m} to obtain the generalized diversity gain 
^2(771) = K^(m) in Theorem [21 
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Figure 1: Two- hop wireless network. 
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Figure 2: Expected rate as a function of number of selected relays. 
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Figure 3: Expected rate for two relay power allocation strategies. 
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Figure 4: Expected rate as a function of average received SNR at destination. 
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Figure 5: Expected rate as a function of power split at relays. 
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